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Abstract 

We apply a self-consistent relativistic mean-field variational "Gaussian functional" 
(or Hartree) approximation to the linear a model with spontaneously and explic- 
itly broken chiral 0(4) symmetry. We set up the self-consistency, or "gap" and the 
Bethe-Salpeter equations. We check and confirm the chiral Ward-Takahashi identities, 
among them the Nambu-Goldstone theorem and the (partial) axial current conserva- 
tion [CAC], both in and away from the chiral limit. With explicit chiral symmetry 
breaking we confirm the Dashen relation for the pion mass and partial CAC. We solve 
numerically the gap and Bethe-Salpeter equations, discuss the solutions' properties and 
the particle content of the theory. 



Address after 1st July 2001: Vinca Institute of Nuclear Sciences, P. O. Box 522, 11001 Belgrade, 
Yugoslavia. 

1 typeset using VTpT^.sty <ver.0.8> 



1. Introduction 



The Gell-Mann-Levy [GML] linear sigma model has long been a subject of nonpertur- 
bative studies, both for its particle physics and statistical mechanics applications. 0*'& In 
this paper we apply a new chirally invariant version of the Lorentz invariant self-consistent 
mean-field variational approximation that goes by many names, inter alia the Gaussian func- 
tional approximation to the linear sigma model. The improvement that we present in 
this paper is the correct implementation of the chiral symmetry in this approximation. We 
prove the chiral Ward-Takahashi identities, among them the Nambu-Goldstone theorem, the 
Dashen relation, and the axial current (partial) conservation (PCAC) in this approximation. 
Then we present a numerically obtained solution of the gap and Bethe-Salpeter equations 
and discuss the particle content of the theory in this approximation. 

Our motivation for this study is the desire to publicize progress made in understanding the 
Gaussian approximation, which is often used in finite temperature/density applications, 
albeit often in incomplete form, and the hope that this work will ultimately lead to the 
clarification of the scalar meson spectroscopy, and in particular of the so-called a meson. 
For this reason, the present study of the GML model must be considered as a methodological 
work providing preparatory for a full-fledged, Nf = 3 calculation. 

This paper consists of six sections. In § 2 we introduce the linear S model. In § 3 we 
outline the Gaussian approximation, and in § 4 we demonstrate its chiral invariance. In § 5 
we present results of the numerical solution of the gap and the Bethe-Salpeter equations and 
analyze the solutions. Finally, we summarize and draw conclusions in § 6. 



§2. The linear E model 

We confine ourselves to the 0{N = 4) symmetric linear a model for the sake of simplicity. 
The Lagrangian density of this theory is 

1 .^2 Trf±2\ 



where 



is a column vector and 



= (00,01,02,03) = (C^,7r) 



n0^) = -l/.^0^+^(0f 

,2 



We assume here that Aq and /ig are not only positive, but such that spontaneous symmetry 
breakdown (SSB) occurs in the mean-field approximation [MFA], to be introduced below. 
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As the chiral symmetry breaking (xSB) term in the Lagrangian, we choose 

^xSB = -Ti-xSB = (2-2) 

as suggested by the underlying NJL quark model. The interaction potential in the new field 
variable s = a — (cr)oB reads 

+ (^%^)(«^ + '^^)^ (2-3) 

The scalar meson (cr) mass, the vacuum expectation value (v.e.v.) and the pion (tt) mass 



are 

s v"^ 

(o-)oB = ^^B = /ttB = 2 + ''^0^ , (2-4a) 

= -/^o + 3Ao/.'b , (2-4b) 
mlB = -lA + Ao/.'b = — ■ (2-4c) 

Note that once the pion mass ttIttb and decay constant /ttb have been fixed, there is only one 
free parameter left in this (tree) approximation, the scalar meson a mass mg-B- The quartic 
coupling constant Aq can be expressed as 

This relation is a powerful result, as it implies a cubic dependence of the cr decay width on 
the its mass: 



r,„ = z^^^^^f^^\ (2.6) 

327r/^B"^<TB V V m„B ) 
Therefore, as soon as the a mass exceeds the two-pion threshold, the decay width increases 
so quickly that instead of a sharp "peak" in the cross section there is a broad "bump", 
unrecognizable as a resonance. This fact is a consequence of the strong coupling implied by 
Eq. (2-5), which in turn is a consequence of the linear realization of the chiral symmetry 
in the Born approximation. Here, natural (and long-standing) questions arise: Does this 
effect survive after taking into account of the loop corrections? What kind of a particle, 
if any, corresponds to the cr field, and how does one identify it? Many studies have been 
devoted to answering these questions, but most suffer from either being perturbative, which 
is unacceptable in the strong coupling case, or from not being chirally symmetric. We satisfy 
both of these requirements by employing the non-perturbative, chirally symmetric method 
described in the next section. 
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53. The Gaussian variational method 



Over the past 20 years we have seen the relativistic Rayleigh-Ritz variational approxima- 
tion based on the Gaussian ground state (vacuum) functional elevated from a little-known 
specialist technical toolll* to a textbook methodic. This method sometimes also goes by 
the names 'self-consistent mean field approximation' (MFA) and 'Hartree + RPA'. In the 
following we use these terms interchangeably. 



3.1. The basics of the Gaussian functional approximation 
We use the Gaussian ground state functional Ansatz 

%[$]=Xexp {-^Jd^ J rfy [0.(x) - (0,(x))]G,r.i(x,y) [0,(y) - (0,(y)))] 



(3-1) 



where A/" is the normalization constant, {(f)i{x)) is the vacuum expectation value (v.e.v.) 
of the i-th spinless field (which henceforth we will assume to be translationally invariant, 

(0,(x)) = (0,(0)) = (0,)), and 

dk 1 



Gy(x,y) = ^6ij J 



^ik-(x-y) 



(27^)' Jk^ + 



Furthermore, note that we have explicitly kept h (while setting the velocity of light c = 1) 
to keep track of quantum corrections and count the number of "loops" in our calculation. 
Then the "vacuum" (ground state) energy density becomes 

£(m„ (0,)) = -e(0o) - l^^l{c|>f + ^ [{cf>)f 

r 1 1 



4< 



where 



hinii) 



dk 



+ 3/1^ ^ j2(^^) + /o(m,)/o(m, 

i i<j 



1 



d'^k 



G'ii(x,x) , 



(2vr) 



(3-2) 



(3-3) 



2^ i^^r.J^+^ (27r)4P-m2 + ze 

1 f dk [— i f d'^k , /.o 2 \ / .\ 

-^k2 + mf = -^J 7^ log [k -mi+tej+ const. (3-4) 



*' Due to the Bose statistics of our fields and the covariance of our approach, RPA might be equivalent 
to the Tamm-Dancoff approximation (TDA) in this case. 
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We identify hli{mi) with the famihar "zero-point" energy density of a free scalar field of 
mass rrii. 

The divergent integrals /o,i(?^i) are understood to be regularized via a UV momentum 
cutoff A. Thus we have introduced a new free parameter into the calculation. This was 
bound to happen in one form or another, since even in the renormalized perturbation theory 
one must introduce a new dimensional quantity (the "renormalization scale/point") at the 
one loop level. We treat this model as an effective theory and thus keep the cutoff without 
renormalization. 



3.2. The gap equations 

We vary the energy density with respect to the field vacuum expectation values {(pi) and 
the "dressed" masses m^. The extremization condition with respect to the field vacuum 
expectation values reads 



0; 



0, 



or 



' dS{mi, {(t>i)y 
' d£{mi, {(t>i)y 



-e + (0o) 



-/i^ + Ao (0)' + 3Mo(mo) + /o( 



1=1 



j^k=0 



= 

. 

(3-5) 



Note that if we assume that (0o) and (0i,2,3) are simultaneously nonzero in the chiral limit 
6^0, then after subtracting one of the equations (3-5) from the other, we are forced to 
conclude that /o(M = mo) = /o(yU = ^1^2,3), or that these two masses are identical. This, 
however, leads one to the symmetric phase of the theory, so we must ignore this possibility. 
Instead, we assume only one to be nonzero, e.g. (0o) = 0, while (01,2,3) = 0. Thus 
the first set of energy minimization equations in the Gaussian variational approximation B 
reads 



e 

— + Ao 

V 



v'^ + 3hIo{M) + Shloii^) 
1,2,3 , 



(3-6a) 
(3-6b) 



where the divergent integral Io{mi) is understood to be regularized via a UV momentum 
cutoff A, either three dimensional or four dimensional. 



*•* There are several renormalization schemes for the Gaussian approximation, but as their names ("pre- 
carious" and "autonomous") suggest, they are unstable. 
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Equations (3-6a) and (3-6b) can be identified with the truncated Schwinger-Dyson (SD) 
equations!^ for the one-point Green function (see Fig. 1). We associate the nonvanishing 




(a) (b) 
Fig. 1. The one-point Green function Schwinger-Dyson equation determining the dynamics of the 
a model in the Hartree approximation: the one-loop graph (a), and the tree tadpole diagram 
(b). The solid line denotes the bare meson multiplet, and the double solid line is the dressed 
meson multiplet. The shaded blob together with the double line leading to it (the "tadpole") 
denotes the vacuum expectation value of the field (i.e. the one-point Green function), and 
the solid dot in the intersection of the four lines denotes the bare four-point coupling. The 
diagrams are explicitly multiplied by their symmetry numbers. 

vacuum expectation value (v.e.v.) with the "sigma meson" field whose apparent mass is 
given by mo = M, and the remaining three fields 0j(2 = 1, 2, 3), of mass mj = /i, form the 
pion triplet. The second set of energy minimization equations reads 



2(0o)' + (0)' + 3hIoiM) + 3Mo(/i) 



= -f^l + Ao [3v^ + 3/1/0 (M) + Shloin) 
= -fil + Ao [{(f>f + hIo{M) + 5Mo(/i) 
= + Ao \v'^ + hloiM) + 5hIo{i2) 



(3-7a) 



(3.7b) 



Equations (3- 7a) and (3- 7b) also have the Feynman-diagrammatic interpretation shown in 
Fig. 2. Upon inserting Eqs. (3-6a) and (3-6b) into Eqs. (3-7a) and (3-7b), the following two 



+ 



(a) (b) 

Fig. 2. Two-point Green function Schwinger-Dyson equation: the one-loop graph (a), and the tree 
tadpole diagram (b). The symbols have the same meaning as in Fig. 1. 




coupled "gap" equations emerge: 



(3-8a) 
(3-8b) 
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One might be tempted to identify jj with the pion mass and M with the cr mass, then solve 
these equations and stop there. However, with e = these equations admit only massive 
solutions M > /i > for real, positive values of Aq and /ig and any real ultraviolet cutoff 
A in the momentum integrals Io{mi) and as these are positive definite (for any real 

mass). In other words the "pion" (0i,02 and (f)^) excitations are massive, with mass /i 7^ 0, 
in MFA, even in the chiral limit. This looks like a breakdown of the 0(4) invariance of this 
method, but, as discussed at length in Ref. |^), there is a simple solution obtained using the 
Bethe-Salpeter equation. Before proceeding to solve the gap equations (3-8a) and (3-8b), 
we will have to determine the value of the e parameter in terms of observables calculated 
in the Gaussian approximation. For this purpose we also have to use the Bethe-Salpeter 
equation. 

3.3. The Bethe-Salpeter equation ("RPA"): an scattering 

In Ref. ^ we have shown that the Nambu-Goldstone particles appear as poles in the 
two-particle propagator i.e. they are bound states of the two distinct massive elementary 
excitations in the theory. We specify the two-body dynamics in the theory in terms of the 
four-point SD equation or, equivalently, of the Bethe-Salpeter equation, see Figs. 3,4. We 



Fig. 3. Four-point Green function Schwinger-Dyson or Bethe-Salpeter equation. The square "box" 
represents the potential, and the round "blob" is the BS amplitude itself. All lines represent 
dressed fields like the double lines in Figs. 1 and 2. 



Fig. 4. The potential (square "box") entering the Bethe-Salpeter equation, as defined in the RPA. 
focus on the s-channel part D.,^{s) of the total four-point scattering amplitude T{s,t,u). Its 



*^ It is well known from the quantum many-body literature that the Hartree or mean-field approxi- 
mation does not respect internal symmetries. The corrective measure goes by the name of random phase 
approximation (RPA). 
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four-point SD equation reads 



(27r)4 [P - M2 + te] [{k - Pf - fi^ + ie] 



with the solution 



s- fi^ 



DJs) 



2Ao 



1 + 



M 



2 £ 



S — 



VJs) 



1 - K(s)i7,(s) 



(3-9) 
(3-10) 

(3-11) 



(3-12) 



where s = {pi +^2)^ = is the center-of-mass (CM) energy. This propagator can also be 
written in the form (see Ref. ITT 



DJs) 



rj2 

s — ml 



(3-13) 



where 



_2 —2 



+ 



2Ao(M2 




M2 



(4vr) 



1 

M2 



M2 



^M2-^2^ 



^M2 -^2^ 



^2Ao(M2 -^2)2 



M2 - /i2 



(3-14) 



We see that the second term in curly brackets is roughly 1% as large as the first one, and 
therefore we may neglect it in the first approximation. 

Here we have simply presented the correct form of the four-point SD equation based on 
the truncation of the exact SD equation. & The derivation from the Gaussian approximation 
in the symmetric phase of the theory can be found in Ref. The corresponding derivation 
in the asymmetric (Nambu-Goldstone) phase can be found in Ref. |^). Furthermore, this BS 
equation is also the "random phase approximation" (RPA) equation of motion that describes 
"quasi-particles" in this theory (see Ref. |6D,P),p!0D). 

3.4. The Bethe-Salpeter equation: tttt scattering 

The dynamics are specified in terms of the SD, or equivalently, Bethe-Salpeter (BS) 
equation for the four-point Green functions Dij{s), where the indices i and j denote the 
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isospin of the pions in the initial and final states, respectively. We use the isospin invariance 
to split this 9x9 matrix equation into three invariant subspaces: (a) isoscalar, (b) isovector, 
and (c) isotensor. Because we consider S-wave scattering, the isovector amplitude vanishes 
identically, due to the Bose-Einstein statistics of the pions. The isotensor BS equation is 
linear and can be solved straightforwardly, but it is without distinguishing features. On the 
other hand, in the isoscalar channel we expect to see the a meson. The corresponding BS 
equations, consist of four coupled equations that can be put into 2x2 matrix form.lHI) 

The main difference between the isoscalar channel and the pion channel, considered in 
§ 3.3, is that here we have two distinct intermediate states, one with two "elementary" sigma 
fields (0o) and another with two "elementary" pion fields = 1,2,3). The isoscalar SD 
equations couple these two channels: 

1 3 

Dmm{s) = Vmm(s) + -Vmm(s)/mm(s)-Da/m(s) + -VA,f^t(s)/^^(s)L)^M(s) , 

DM^l{s) = VMf,{s) + ^VMf^{s)I^,f,{s)D^^{s) + ^VMM{s)lMMis)DMt,{s) , 
DfiM^s) = V^m{s) + ^Vf,M{,s)lMM{s)DMM{s) + ^V^^{s) I ^^{s) D f,M (s) , 

Df,^,{s) = V^^{s) + ^Vf,f,{s)I^f,{s)D^^{s) + ^V^m{s)Imm{s)Dm^{s) . (3-15) 
The equations in (3-15) can be cast into matrix form 

= V + ivnD,, , (3-16) 



Dmm Dm^i 



i D 

3 -^MPi 



(3-17) 



where the boldfaced symbols are matrices. The solution to the matrix equation (3-16) is 



where 



= n - n)-i V 



V 



Vmm 

1 

M 3 

/ 



2An 



1 + 3 



Af2 



-M2 



V 



1 + 3 



M2- 



-M2 



1 + 3 



5 + 3 



s-Af2 



y_ 



(3-18) 



(3-19) 
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and 



n 



'MM 

3/ 



(3-20) 



The invariant functions Iii{s), Jii{s) are given by 



. f cTI 
J 72^ 



(27r)4 [A;2 - mf + ie] [(k - - mf + ie] 



lu{0) + 



{Any 



[-1 + Jii{s)) , 



(3-21) 



where s — P^, and the real and imaginary parts of Jm(s) are 



I Ami 



— 1 arcsin 



4mf 



4mi 



log 



4m,- 



+ 



4m,- 



I — 

2 



(0< 
(1< 



4mf 

s 



<1) 



< oo) 



(3-22) 



The equal mass integrals at zero momentum /^(O) are cutoff dependent constants ("sub- 
traction constants" in the dispersion relations language) whose values are determined by the 
solutions M,fi, and A to the gap equation (5-1). Thus, fixing of the subtraction constants is 
one of the primary consequences of the self-consistent gap equations. We use 



if)(0) = 2(47r)-^ 



X4 



- l0g(x4 + 1) 



Xi 



m 



- ln(v^ + Vl + X3) 



; X3 



(3-23a) 



A \2 

— ) . (3-23b) 

m / 



/X3{1 + X^) 

Thus self-consistency of the gap equations enters into the scattering problem and constrains 
the remaining free parameters. 

We see from Eq. (3-15) that in order to calculate the Dmm amphtude, we need to know 
the Dfj_M{s) amplitude, and vice versa. In other words, we must solve the system of four 
coupled equations (3-15). It turns out that this system splits into two systems with two 
unknowns, with the same discriminant T>. This fact, (a) ensures that there are at most two 
poles in the solutions, and (b) greatly simphfies the algebra. The solutions to the equations 
(3-15) are 

1 



Dmm{s) 



V{s) 
1 



V{s) 



(Vmm(s) - 12Ao/^^(s)V^mm(s)) , 



Vmh{s) 



V{s) 



(3-24) 
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where 

V{s) = 1 - ^ [Vmm{s)Imm{s) + V^^{s)I^^{s)] + 6\olMMis)I^^{s)V^M{s) (3-25) 

is the discriminant of "one half of the system of equations (3-15). Elementary and composite 
states in the s-channel manifest themselves as poles in the Da-(s) matrix, or equivalently as 
roots of 

{s - M^)V{s) = 0. (3-26) 

This equation is identical to the formula for the a mass in the optimized perturbation theory 
(OPT)© 

An inspection of Eq. (3-25) leads one to think that (s — M^)'D(s) should be a quadratic 
polynomial (function) of s and therefore have two roots, since each {s-M^)Vij{s) is a linear 
function of s. However, as a consequence of chiral symmetry, one finds that the function 
(s — M^)P(s) is (only) linear in s. This, in turn, ensures that there is only one root of the 
"eigenvalue" equation (3-26), and hence only one pole in the a propagator at least for small 
values of s, where the la functions' logarithmic s dependence may be neglected. For large 
values of s, this is no longer the case, as we show in § ||. 



§4. Chiral symmetry Ward identities in the Gaussian approximation 

Chiral Ward-Takahashi identities follow from the underlying chiral symmetry of the linear 
sigma model and typically relate (n-l)-point Green functions to other n-point functions 
and/or currents. These identities were developed by Lee at the perturbative one loop level, 
i) and by Symanzik at arbitrary orders of perturbation theory, We call them the Lee- 
Symanzik [LS] identities. To our knowledge, for selected infinite classes of diagrams in the 
linear sigma model no proofs of LS identities were given prior to Ref. |^). The NG theorem 
is the simplest LS identity. When the chiral symmetry is explicitly broken, NG theorem 
turns into a relation between the chiral symmetry-breaking parameter and the NG boson 
mass, as first discussed by Dashen.!!!* The NG theorem in the chiral limit has already been 
addressed in the Gaussian approximation and related formalisms in Refs. ^,|D and |16D. For 
this reason, we consider the nonchiral case. 

4.1. D ashen' s formula and the Namhu-Goldstone theorem 

As shown in Ref. |^) in the chiral limit the Nambu-Goldstone particle appears as a zero- 
mass pole in the pion channel two-particle propagator "Ec^. (3-12). Next, we consider the zero 
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CM energy P = polarization function V.,r{0)n.^{0) in the nonchiral case. We use Eq. (3-7b) 
to write 



v;(o)i7.(o) 



2\nh 



(M2 - /i2) 



and then use Eq. (3-7a) to obtain the final result, 



[Jo(M) - /o(/i)] 
M2 - ■ 

1 TT- 



K(0)i7.(0) 



e 

V 



M2 



The propagator Eq. (3-12) evaluated at zero momentum can be written as 



(4-1) 



(4-2) 



^.(0) 



„2 



-mt 



2Ao' 
ml 



M2 



v;(o) 



M2 



M2 



M2 



(4-3) 



1 - V;(0)i7,(0) 

which leads to the result e = m'^v + 0{e'^). This is also in agreement with a general result due 
to Dashen.0'' The vacuum expectation value of the U operator, when sandwiched between 
two vacuum states, yields Dashen's formula. 



{fm'fY' = famlJ, = -(0 I [QI [Ql H^sb]] \ 0) + 



(4-4) 



for the pseudoscalar meson mass squared {rripg) and decay constant fa for an arbitrary chiral 
symmetry-breaking term in the Hamiltonian density Ti^SB- This is a model-independent 
result based on the equations of motion in the Heisenberg representation and on the LSZ 
reduction formulas. We use the canonical commutation relations and the axial charge to 
evaluate the E operator in the linear E model. 



E5' 



ab 



QUQi'>ixSBm 



-ea5' 



ab 



(4-5) 



which vanishes when the chiral symmetry is not explicitly broken. Taking the vacuum 
expectation value of this expression, we find 



(4-6) 



which leads to 

e = mlU . (4-7) 

This relation is satisfied both in the Born and Gaussian approximations, which indicates 
that both are chirally symmetric. Now that e has been fixed, note that the gap equation 
(3-8b) implies < //^^ in agreement with the variational nature of the Gaussian approxi- 
mation. This is important in the numerical calculations discussed below and for the physical 
interpretation. 
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4.2. Axial current Ward identity 



n 




Fig. 5. Axial current matrix element 

The axial current matrix element corresponding to the Feynman diagrams displayed in 
Fig. ^ reads 

JUip\p) = mp')\j,mMp)) 

= (P + P)p. + % I o _ I 



— /i^ 



^M5(g)^7r(g) 



where 1)^5(0') is defined by 



(27r) 



{2k + g)^ + 



'M 



2 £ ■ 



— /i^ 



- M2] \{k + g)2 - /i2 



(4-8) 



(4-9) 



and satisfies the chiral Ward identity. 

,2 



18D 



2X0 



(K(0)i7^(0) - K(g')i7.(g')) - - (i7.(o) - n^q 



(4-10) 



We insert the vertex r^^{q) in Eq. (4-9) together with the two-body propagator DT,{q^) in 
Eq. (3-12) into Eq. (4-8) to find 

JAp ^P) = (P +P)f^ + % 



(f — ^ 



r,2 



X 



^ (K(0)i7.(0) - V;(g')i7^(g2)) - - (i7,(0) - i7^(g- 



Vl-V;(g2)i7,(g2); 



g2 _ ^2 



(4-11) 



Here we have q'^ = (p — pY . This result manifestly lacks a pole at g^ = /i^. The composite 
state plays precisely the role of the Nambu-Goldstone boson in the conservation of the (axial) 
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Noether current, i.e., in the basic axial Ward-Takahashi identity 

cfJ..{p,p) = (p'^ - - - M^) + ml - l) , (4-12) 

which follows directly from Eq. (4-11). Furthermore, /,r is defined by 

(0| J^5|i7(q)) = Uq)% = 9'^^r,,iq), (4-13) 
from which (in the chiral limit) it follows that 

MO)g^" = M' - (4-14) 

by way of the axial Ward identity Eq. (4-10). This result, together with Eq. (3-14) for g^^^ 
forms the basic result for the composite pion decay constant, 

fn - MO) = g-Jf (m^ -ii'')=v[l + ((47r)-2)] V. (4-15) 

§5. Numerical solutions 

5.1. The self- consistency or gap equation 

Having determined the value of the parameter e — vml in terms of observables, we can 
solve the gap equations. We fix the v.e.v. at a value of 93 MeV for the pion decay constant 
/^r and a value of 140 MeV for the physical pion mass m^. As a result, the system of gap 
equations (3-8a) and (3-8b) turns into a single equation: 

v' = fl = i^^]n[iM - hm 



III 



2 ' 



^2 _ 



^ h[lM-Io{M)] . (5-1) 



Here, we have used Eq. (3-3) as the integral to be regulated. We give here the results for 
(covariant) four-dimensional Euclidean cutoff regularization the three-dimensional regular- 
ization of this quadratically divergent integral: 



/^(m') = (47r)-W [X4 - ln(l + xA] ; X4 = f — V , (5-2a) 
r(3)^_2, _o^.^^-2_2 .fZJV^::r,_^„f /:;r^./T^^^. xs^(^)^ . (5-2b) 



^Jx3{l + X3) - ln(v^ + Vl + X3) 



\ m . 

Note that the ^ — > limit of Eq. (5-1) is non-trivial: on the right-hand side both the 
numerator and the denominator vanish. Further, note that one parameter, in this case the 
cutoff A, remains free. 
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Numerical solutions to the equation (5-1) are plotted in Figs. 6 and 7 for various values of 
the cutoff A. Every point on the (/i, M) curve represents a solution to the gap equation, thus 
signalling the existence of freedom of choice in the form of one continuous free parameter. 
This free parameter can be related to the bare coupling constant Aq by Eq. (3-8a) for every 
(/i, M) pair. 

In Figs. 6 and 7 we see that as the cutoff A increases, all solutions to the gap equation 
approach the symmetry restoration limit M fi for large values of M, or equivalently large 
values of Aq. This implies that the large boson loop effects lead to symmetry restoration, in 
contrast to the fermion loops, which lead to symmetry breaking. Solutions that lie above 
the 2/x threshold require rather small values of (either kind of) cutoff A. However, this is in 
agreement with the small (second) meson-loop cutoff found in l/N^ studies of the NJL chiral 
quark model. 0^ This does not mean that the loop effects are necessarily small, however. 



5.2. The Bethe-Salpeter or scattering Equation 

We reduced the BS equation (3-15) in the isoscalar channel to solving a single alge- 
braic equation (3-26) involving the transcendental analytic functions Imm{s) and I^^{s) with 
branch cuts at imaginary parts above the corresponding thresholds. This equation has, 
in general, real and imaginary parts: For a mass values lying below the two-body threshold, 
only the real part is relevant, while for heavier a masses the imaginary part must be taken 
into account as well. The latter determines the natural decay width of the a meson. 

From numerical solutions to the real part of Eq. (3-26) shown in Fig. 8, it can be seen 
that the a mass is always shifted downward from the "elementary" sigma field (0o) mass 
M, in agreement with the variational property of the mean-field approximation. For small 
values of M (< 2yu) (i.e., when the coupling constant Aq lies below some critical value Ac, 
which is a function of the masses fi, M and the cutoff A), the a meson mass (i.e., the real 
part of the pole position) drops below the 2/i threshold, and the a meson comes to consist 
predominantly of the bare 0o state with some 27r and 2a "cloud" components admixed to 



its wave function, With increasing coupling constant Aq, the physical a's mass increases 
above the 2fi threshold, and the "bare" and "dressed" components of the wave function 
can no longer be separated. Then, the state itself must be considered as predominantly 
a meson-meson composite. For weak couplings, only one state has been found to exist in 
the a channel of the MFA. In Fig. 8 it can be seen that the a mass changes continuously 



*^ Imm{s) and I^^{s) have logarithmic branch points and therefore infinitely many sheets, in contrast 
with the nonrelativistic case, in which the branch points are of the square root type, with only two sheets. 

**^ Another interpretation of these results (that might be only seniantically different from that one) has 
been given in the language of operator many-body ( "quasi-particle RPA") methods. H^'EH) 
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with decreasing coupling Aq and connects smoothly to the perturbative a mass in the weak 
coupling limit. 




Fig. 8. The solution rua to the Bethe-Salpeter equation in the isoscalar channel as a function of 
the variational parameter M for various values of the cutoff, A. The curves denoted 2/i show 
the movement of the 2/x threshold for corresponding values of the parameters M and A. 

Note, further, that for many values of the cutoff A and above some critical value of M, 
far into the 2/i continuum, there is a second root of the real part of Eq. (3-26). As increases 
M ~ a/Aq, the two roots sometimes merge into one and then immediately disappear (see 
Fig 9). For other values of the parameters, the two roots diverge, the smaller one moving 
down to zero, while the heavier one moves back up in mass again. In either case, the smaller 
zero, which is connected to the perturbative solution, has an upper limit generally below 1 
GeV. This is perhaps the most interesting result of this paper. We remind the reader that 
these results are not artifacts of the Gaussian approximation viewpoint, as exactly the same 
equations govern the OPT meson masses. 

The question of the physical meaning of the two zeros arises; i.e., if actual poles exist 
on the second unphysical Riemann sheet of the S'-matrix, that can be associated with these 
zeros in the real part of the inverse propagator? The larger zero is almost certainly not a 
conventional pole, because the derivative of Eq. (3-26) evaluated at the root has sign opposite 
to that at the smaller root. Because this derivative is related to the effective couphng constant 
squared, this implies that the upper state has a non-Hermitian coupling to the bare states. 
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Fig. 9. The solution rrifj to tlic Bethe-Salpctcr equation in the isoscalar channel as a function of 
the variational parameter M for a small cutoff, A = 0.2 GeV. Note the double-valued nature 
of this functions, as well as the absence of solutions above certain values of M. 

This question is more difficult to address, as it demands analytic continuation onto the 
second Riemann sheet. It will be left to a future investigation. 



§6. Summary and conclusions 

In summary, we have: 1) constructed a unitary, Lorcntz and chirally invariant, self- 
consistent variational approximation to the linear a model; 2) solved the coupled self- 
consistent equations of motion in this, the mean-field plus random-phase approximation 
(MFA + RPA). The solutions to these equations also determine the minimum of the "op- 
timized perturbation theory" effective potential.; 3) shown that the particle content of the 
mean-field plus random-phase approximation to the 0(4) linear sigma model is the same as 
in the Born approximation, at least for weak coupling, i.e., there are three Goldstone bosons 
TT and one a state. 4) found that the pions' mass is unchanged, as a consequence of the 
vahdity of Dashen's relation in the MFA -|- RPA, whereas the a mass and width can differ 
substantially from the Born values, depending on the free parameters. 5) calculated the a 
meson mass by solving the nonperturbative Bethe-Salpeter equation. We found a second 
solution for large values of mass and coupling, whose physical interpretation is yet unclear. 
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The mean-field or Gaussian method was initially fraught with problems when applied to 
the linear a model with spontaneously broken internal symmetry - the Goldstone theorem 
did not seem to "work." This problem was solved in Ref. |^: The Goldstone boson found in 
the Gaussian approximation & turns out to be a composite massless state, just as in the NJL 
model. Yet, there seemed to exist another massive state with the quantum numbers of the 
pion. In fact, however, this is true only in appearance: There is no pole in the propagator 
corresponding to this "particle". The MFA to the bosonic linear a model is significantly 
different from the NJL one in one regard: Whereas in the NJL model, the gap equation 
describes "dressing" of the fermions and the BS equation describes mesons as bound states 
of dressed fermions, in the linear a model both the gap and the BS equations describe 
(two different) "dressings" (first and second renormalizations) of mesons. Therefore, the 
results of the intermediate (first) renormalization remain in the theory even after the second 
renormalization and confuse the issues of the physical content. A more complicated situation 
exists in the scalar sector. We have not yet written the last word on this subject. 

We hope to extend these calculations to physical applications for more realistic La- 
grangians in the future. 
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